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Abstract 

We argue that string percolation is in the origin of i) an approximately flat 
rapidity distribution and of ii) an approximately constant forward-backward 
correlation parameter b over a substantial fraction of the available rapidity 
Predictions are given for pp collisions at LHC, y/s = 14TeV and y/s = 
5.5TeV. 
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In AA (nucleus-nucleus or hadron-hadron) high energy collisions one nat- 
urally expects positive Forward-Backward (F-B) correlations to occur. This 
happens because the overlap area of interaction acts both in the forward and 
in the backward direction. If the impact parameter is small, many parti- 
cles are emitted forward, and many particles are emitted backward. On the 
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contrary, if the impact parameter is large, few particles are emitted forward 
and few particles are emitted backward. Note that in addition to these F-B 
correlations, local correlations in rapidity exist due, for instance, to reso- 
nance decays. As these short-range correlations are not supposed to depend 
on impact parameter, it is clear that F-B correlations increasingly dominate 
with the decrease of impact parameter, or the increase of the number of 
participants N part _. (See, for instance [1]). 

The observation of F-B correlations is then not a surprise. What was 
not anticipated was that the strength of the F-B correlations was practically 
constant over a large rapidity interval, perhaps of the order of the beam 
rapidity interval AY. Recent results at RHIC (including the discovery of the 
ridge phenomenon) show that the F-B correlations extend over a large region 
in rapidity [2,3,4,5]. 

Monte Carlo event generators, such as HIJING, predict a fast decrease of 
the F-B correlation with the rapidity (pseudorapidity) interval At] [6]. The 
Parton String Model [7], which includes pair of strings fusion, shows a slower 
drecrease of the correlation with Arj, thus being in closer agreement with 
data (see [8]). 

Models that can explain, in a straight forward manner, long range F-B 
correlations are models that introduce extended objects in rapidity as colour 
flux tubes [9] or strings [10]. 

In conventional F-B correlation studies b is the F-B correlation parameter 
defined by the relation [11] 

(n B ) (n F ) ~ bn F , (1) 

where n F is the number of particles in the forward rapidity window and (n B ) 
the average number of particles in the backward window. The correlation 
parameter b is given by 

b = [(n F n B ) - (n F )(n B )\ / [(n 2 F ) - (n F ) 2 ] . (2) 

Adopting the two step scenario [12,13], where extended objects are formed 
first, followed by local emissions of particles, assumed Poissonian, one obtains 

, _ (n B )/(n F ) 

b -l + K/(n F )> {3} 

Where K is the inverse of the normalized variance of the number of collisions 
v distribution: 
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For symmetric windows (n B ) = (rip) and (3) becomes 



It is now clear that in this formalism an almost Arj independent F-B 
correlation means that the particle density distribution is approximately flat 
around mid rapidity. Particle densities and F-B correlations are a kind of di- 
rect image of the extended objects stretched between the interacting hadronic 
sheets [9]. 

In the dual string model (see, for instance, [1]) strings are constructed 
from partons, one from each nucleus or hadron, with Feynman x values x + 
and respectively. At high energy or density most of the partons are 
wee partons (gluons) with, in general, x + ~ x_ ~ 1/y/s, where -/s is the 
centre of mass energy. In such situation the formed strings are small, with 
fixed length Ayi, in rapidity. One does expect the particle distribution to 
be peaked around mid rapidity, and being not flat at all in a large rapidity 
interval. But the region of x + ~ x_ ~ is the region of saturation and 
percolation of strings. If (N s ) strings percolate the size in rapidity of the 
percolated strings is [14,15], 

Ay {Ns) = A yi + 2\n(N s ) . (6) 

In the naive model with just these percolating strings, if we impose energy 
conservation [15] we obtain 

(7) 

with A = 2/7, such that 

Ay {Ns) = 2XAY ~ l -AY . (8) 

This simple model tells us that half of the full rapidity interval AY is occupied 
by a flat distribution. 

One can be more sophisticated and use an improved version of the step 
function, the Fermi-like distribution [16, 17] that fits RHIC/PHOBOS data 
and can be obtained by an evolution equation for the particle density [18], 

2 dn e XY 
-■■ B 7' (V>0), (9) 



Npart. dy + i 
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where N part is the number of participating nucleons, Y is the beam rapidity 
and A, a and 5 are parameters (from fit to RHIC data: A = 0.26, a = 0.31 
and 5 = 0.75). The region of approximately flat particle density is between 
the inflexion point, at i] > 0, and the corresponding point at i] < 0: 

Ar) Flat ~ (1 - a)AY ~ 0.7AF . (10) 

Eq. (9), however, strictly speaking, is only valid for large values of N part . 
Having in mind, on one hand, the presence of nucleon-nucleon contribution 
and, on the other, multiple nucleon contributions, we write, in the high 
energy limit, see [19], 

2 dn l._, .2 .1.. e XY 

~ o( (-*- ~ ("at ) 5 )) ,-(i- a)y > (11) 



Npart. dy 2 ^Vport. e~ *^ + 1 

For large values of iVp art . we recover (9). For mid rapidity the validity of the 
low energy version of (11) was tested in [20]. 

Recently, in [19], the two step scenario was successfully applied to central 
Au-Au collisions, RHIC data at y/s = 200GeV [21]. In particular the correla- 
tion parameter b was predicted to be fairly constant up to Ar/<3. Predictions 
were also presented for Pb — Pb at LHC, with an approximately constant b 
extending over a region A?7<6. 

Recent data on correlations for pp collisions at RHIC [5] have shown that 
the rapidity correlation is peaked at i] = - which is typical of a non- 
percolating regime - and that no azimutal elliptic flow effects, contrary to 
central Au-Au, were detected. However there are strong indications that the 
physics of pp collisions at very high energy is the same as the physics of heavy 
ion collisions at lower energy [22], and we then expect that the formalism of 
[19] can be applied to pp collisions at LHC, y/s = 14TeV. 

In fact, in percolation theory the transition to the percolation regime 
occurs for a critical value of the transverse density, 

Vaa=(^) 2 (N s ) (12) 

where r is the single string radius, and R, the effective radius of the overlap 
region of interaction. Simple geometrical or Glauber like arguments give R = 
Rp(^Y^)^ and (N s ) ~ (iVf)(^^)t, where (N s ) is the average number of 
strings and (iVf) the average number of strings in pp collisions. In [10] it was 
estimated that in central PbPb collisions the transitions with rf PhPh — 1.15 
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Figure 1: The F-B correlation parameter b at LHC, y/s = 5.5TeV, in PbPb central colli- 
sions. 

occurred at SPS for y/s ~ 20GeV. By assuming that (Nf) ~ sr we estimate 
rjp p to occur at LHC energies, for \/s pp ~ 6TeV. 

In order to estimate b, see (5), one needs (rip) and K. The window 
multiplicity is estimated by using (11), under the assumption of a window 
with width 5rj = 0.2 (as in [19]). The problem is to estimate K, (4), (in [19] 
K was a free parameter). 

In a simple model - again a two step scenario - where nucleon-nucleon col- 
lisions are followed, in an independent manner, by parton-parton collisions, 
we can write [23], 



where v 9 ( z/ ) is the probability of having v collisions, ^{vn) the probability 
of having nucleon-nucleon collisions and p(/i) the probability of having, 
in a nucleon-nucleon collision (say, pp), /i partonic collisions. From (13) we 
conclude, 



(13) 



i=i 



(u) = (u N )(fj) , 



(14) 



and 
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Figure 2: The F-B correlation parameter b at LHC, yfs — 14TeV, full line, and y/s = 
5.5TeV, dashed line, in pp collisions. 

(^ 2 ) - H 2 H) - (^v) 2 , (^ 2 ) - (/^) 2 i n ^ 
<^> 2 <^> 2 (^) 2 M ' 1 } 

or, 

1 111 

where ^ = (fi) 2 /((fi 2 ) - (^u) 2 ), and K N = (v N ) 2 /{{v 2 N ) - ("n) 2 )- It is clear, 
from (16), that K (for the overall collision distribution) is smaller than K N 
(for the nucleon collision distribution). For N — 1 (pp collisions) X = fc^. 
As (rip) ~ dn/dy ~ Afp ° rt ' , (9), the correlation parameter 6 increases with 
the N part , we have, in addition to (16), the relation (see (5)), 

A < V (17) 



Npart 



Note that the quantities -ft^v and (v N ) are fixed by Glauber calculus [24]: 
Xjv ^ 300 and (u N ) ~ 800. 

By using (16) and (17) we are now in conditions of finding the parameters 
at LHC energies: K ~ 170 for central PbPb and ~ 0.5 for pp collisions. 
The results for b as a function of Ar] are shown in Figs. 1 and 2. In Fig. 3 
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Figure 3: Particle density, dn/dr] in pp collisions at yfs = 14TcV. Errors due to a 3% error 
atributed to the parameter A. 

we present our prediction for dn/dy for pp collisions at \fs = lATeV. In Fig. 
4, we show our expected dependence of 2/N Part x dn/dr] on N ParL (11) at 
LHC, v/i = 5.5TeV. 

We would like to finish with two remarks: i) Even in pp collisions is 
possible to test the impact parameter dependence of the correlation param- 
eter b, by comparing events with different multiplicities (bellow the average 
and above the average, for instance). We expect higher multiplicity events 
to correspond to higher values for b. ii) In [25] it was made a prediction for 
dn/dr](pp, y/s = 14TeV)70±8 charged particles. We expect (see Fig. 3) 80 ±7 
slightly larger than the value obtained in [25]. 

The possibility of existence of collective effects in pp collisions at very 
high energy was considered in [26]. 
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Figure 4: 2/Np art . x dn/dr] versus Np art . at LHC, y 7 ! = 5.5TeV. 
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